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We reply to the criticism from the authors of arXiv: 1502.02782 of the spin-fluctuation scenario 
for charge order in the cuprates. The authors of of arXiv: 1502.02782 argued that spin-fluctuation 
exchange cannot give rise to charge order with observed momentum (Q, 0)/(0, Q ) due to the absence 
of nesting for a half of fermions involved. We explicitly show the instability towards charge order 
exists even in the “worst-case” scenario of anti-nesting for a half of fermions. 


In the recent preprint [j] the authors performed the 
numerical study of superconducting and charge-density- 
wave instabilities in the model of fermions interacting 
via magnetic fluctuations. They argued that there is no 
instability in the charge-density-wave channel with mo¬ 
mentum (0, Q)/( 0, Q) (CDW-x in their notation) because 
there is no analog of Cooper logarithm due to the lack 
of Fermi surface (FS) nesting for one half of the pairs 
of hot spots. This result contradicts earlier works [2j-l4|], 
including our earlier analytical study [4]. In this reply, 
we consider the same model as the authors of 0 and, 
like them, calculate and compare the eigenvalues in the 
superconducting (SC) and CDW-x channels, A SC (T) and 
A c dw(T) (defined such that the system develops an insta¬ 
bility once A(T) becomes greater than one). To make our 
point, we only consider the “worst-case” scenario of per¬ 
fect anti-nesting for one pair of hot spots. We show ex¬ 
plicitly that A c dw{T) diverges logarithmically at T —>• 0, 
in the same way as A sc , and the ratio X c d w /X sc is just a 
constant (smaller than one), whose magnitude is deter¬ 
mined by the ratio of spin-fermion coupling and Landau 
damping. As the consequence, A c dw crosses A c dw = 1 at 
some finite T c d w . We believe that Ref. [0 missed this 
instability because the strength of the Landau damping 
was taken to be much larger than one would get by cal¬ 
culating it self-consistently. 

We also comment on the validity of neglecting the 
fermionic self-energy, as it was done in Ref. [ 0 . We ar¬ 
gue that keeping self-energy is essential even in the non¬ 
quantum-critical, Fermi liquid regime for proper treat¬ 
ment of the effects due to FS curvature. We also com¬ 
ment on the claim in [0 that other potential difference 
between their and earlier works is that previous analy¬ 
sis relied on an expansion around the hot spots, while 
the calculations in Ref. jl] were performed by integrat¬ 
ing over the full Brillouin zone. We argue that for the 
parametrs used in [I], the contribution to the momen¬ 
tum integral comes predominantly from the vicinity of 
hot spots, and thus the expansion around hot spots is a 
good approximation. 

To make explicit comparison with Ref. [ 0 , we con¬ 


sider the same model as they did, namely fermions in¬ 
teracting by exchanging magnetic fluctuations peaked at 
K = (7r,7r), and, like them, we neglect fermionic self¬ 
energy. The dispersion considered in Ref. T] is repro¬ 
duced in Fig. 1. For this dispersion, the velocities at hot 
spots separated along X or Y direction are either almost 
anti-parallel (for one half of hot spots) or almost paral¬ 
lel (for the other half). For SC channel, the velocities of 
relevant fermions are always anti-parallel. The form of 
spin-mediated dynamical interaction used in Ref. [lj is 
the same as in earlier works 30 : 
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where g is the effective spin-fermion coupling constant, 7 
is the Landau damping coefficient, and £ is the magnetic 
correlation length. 

To simplify the presentation, we take velocities at hot 
spots to be either strictly parallel or antiparallel, neglect 
FS curvature at hot spots, and compare the eigenvalues 
in SC and CDW-x channels by analyzing ladder equa¬ 
tions for the corresponding vertices (see Fig. 2). The 
SU(2) symmetry of the hot spot model (Refs.jbf, i7|) im¬ 
plies that the kernel in the SC channel and the one in 
CDW-x channel for anti-parallel velocities are identical. 
The difference between two channels is in that for CDW- 
x channel the interaction peaked at (7 r, 7r) takes fermions 
from the region where Fermi velocities are anti-parallel 
(region 1-2 in Fig. 1) to the region where Fermi veloc¬ 
ities are parallel (region 3-4 in Fig. 1). As the conse¬ 
quence, one needs a composite second-order process to 
bring fermions back into the region where Fermi veloci¬ 
ties are antiparallel, see Fig. 2c (a similar reasoning has 
been recently applied to the analysis of the optical con¬ 
ductivity in the cuprates 0 ). 

The integration over two Green’s functions with an¬ 
tiparallel velocities gives rise to a conventional Cooper 
logarithm log (w s f/T), where w s / ~ £ - 2 / 7 , and the loga¬ 
rithm comes from the smallest momenta and frequencies. 
To logarithmic accuracy, the eigenvalue/eigenfunction 
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equations for superconducting and CDW condensates, < b sc (U c ) and & c dw(k x ) are (small k x is along the FS near 

hot spots 1 and 2 ) 

I 


\sc{T)<f>sc{k x ) = 3^1^ log (u s f/T) f dk' x $ sc (k' x )x(k x - k' x , 0) 

A cdw{T)$ cdw (k x ) = log {usf/T) f dk' x $ cdw (k' x )xcom(k x , k' x , 0) (2) 


where the numerical prefactors are the results of summa¬ 
tion over spin indices and the shift by K is absorbed 
into the shift of momenta. The interaction is attrac¬ 
tive in both SC and CDW channels if we additionally 
assume that both 4> sc (k) and 4) cdw {k) change sign under 
k —>■ k + K (Refs. @0). 

The composite Xcom{k x , k x , 0) is the convolution of two 
dynamical spin susceptibilities and two Green functions 
of fermions with parallel velocities: 


( 3 ) 

X(kx ~ Px,Py,Cl)x(K ~Px,Py,Cl). 


Xcom (k x , k x , 0 ) — 

3 g f dVlmdpxdpy 
87 r 3 J (iCl m - v FPx y 

where momenta p x , p y and frequency Cl are for fermions 
in region 3-4 where velocities of the two hot fermions are 
parallel. If we approximated y(q, Cl) by its static part, 
or took it on the FS (i.e., set p x = 0 in the susceptibil¬ 
ities), the integral would vanish because of double pole 
in l/(iCl m — vfPx) 2 ■ This is what the authors of Ref. [I] 
probably meant when they argued that there is no CDW- 
x instability due to anti-nesting for a half of hot spots. 
However, if we use the full forms of %(q, Cl), we immedi¬ 
ately obtain that the integral in <[3|) does not vanish and 
yields [f| 


Xcom(k x ,k’ x ,0) = Ax{k x - k x ,0)f(k x £,k x £), A = 
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where /( 0 , 0 ) = 1 and to a good numerical accuracy, 
f{x,y) remains close to one for relevant x,y = O(l). As 
a result, Xcom{k x ,k' x , 0) and x(k x - k x , 0) differ just by 
a constant A. Using we then find that the ratio 
Xcdw/^sc = A is independent on £. This result implies 
that A cdw has the same logarithmic dependence on T as 
A sc and therefore must cross one at some finite T cdw . 
This temperature is numerically, but not parametrically 
smaller than superconducting T c . 

The magnitude of A depends on the interplay between 
spin-fermion coupling and Landau damping coefficient. 
In self-consistent calculations within the spin-fermion 
model 7 by itself is expressed via g as 7 = 4g/(irVp). 
(Refs. [5j, llj), hence A = 3/(8n) is just a number. The 
value of A increases if one additionally assumes that 
CDW-x order emerges from some pre-existing pseudo¬ 


gap state which reduces Landau damping 101-112] due to 


a reduction of a low-energy fermionic spectral weight in 
the hot regions. The Landau damping coefficient used in 



FIG. 1. The Brillouin zone, the magnetic Brillouin zone, and 
FS using the same dispersion as in Ref. [J. Hot spots 1,2 
and 3,4 are defined as points on the Fermi surface that inter¬ 
sects with the magnetic Brillouin zone. The red wavy lines 
represents the interaction mediated by spin fluctuations of 
momentum (n,n). 


Ref. [Ij was not obtained self-consistently and was taken 
to be 20eV _1 for g = 0.9eV and v F = 0.78eV, which is 
about ten times larger than self-consistent value. This 
yields a much smaller A = 0.01, which is probably the 
reason why CDW-x instability has not been observed in 
Ref. 0. 

We next discuss the effects of fermionic self-energy and 
FS curvature. The self-energy was neglected in Ref. 
[I] (and in this respect their approach is not a truly 
strong coupling). It plays the central role in relating SC 
and CDW instabilities in the quantum-critical, non-FL 
regime when relevant T > co s f- This regime requires a 
separate consideration Q and we will not discuss it here. 
At energies below u> s f, which we consider here, the self¬ 
energy has a Fermi liquid form and gives rise to mass 
renormalization m*/m = 1 + A, where A = 3g^/(4 ttv f ). 
The self-energy is not overly important for Eq. © as 
it only adds 1/(1 + A) to the overall factor in the r.h.s. 
but does not affect the logarithm. The inclusion of the 
self-energy is, however, essential for the analysis of the 
distractive effect of FS curvature on T cdw . One can 
easily check that the curvature (the term nk x in the 
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FIG. 2. The superconducting vertex <3> sc [Panel (a)] and the 
CDW-x vertex Q c dw [Panels (b) and (c)[ in the ladder approx¬ 
imation. Panel (b) shows the relation between CDW conden¬ 
sate formed by fermions at hot spots 1 and 2 and the one 
formed by fermions at hot spots 3 and 4. In Panel (c) we 
integrated out fermions near hot spots 3 and 4 and obtained 
the self-consistent equation for the CDW condensate formed 
by fermions at hot spots 1 and 2. 


fermionic dispersion near hot spots 1 and 2 ), cuts the 
Cooper logarithm in the CDW channel at a frequency 
of order re£ _2 /(l + A). Using k ~ Vp/Ep, we find that 
this frequency is of order u s f(g/Ep)/( 1 + A). The loga- 
rithmical dependence on T in the CDW channel survives 
as long as (g/Ep )/{[1 + A) 1. Without 1 + A term 
coming from self-energy, this condition is satisfied only if 
g Ep, while for parameters in Ref. 0 g ~ Ep. As a 
result, without self-energy, the FS curvature strongly re¬ 
duces T c d w . With self-energy, the logw s //T behavior in 
CDW-x channel definitely holds at large enough £, when 
A 1. As the consequence, at large enough £, T c d w is 
only weakly affected by the FS curvature. As £ decreases, 
the effect of the curvature gets stronger, T c d w decreases 
and eventually vanishes at some finite £, as we argued be¬ 
fore 0] The same behavior due to FS curvature holds for 
charge bond order with momentum ( Q,Q ) (Refs. 00 )- 
Finally, we discuss the validity of the expansion around 
hot spots. From Eq. (4) we see that the relevant momenta 
k x , k' x (which, we remind, are deviation from hot spots 1 , 


2 along the FS) are of order £ _1 . The typical momenta 
p x and p y in Eq. (3) are also of order £ . Taking, e.g., 

£ = 3a, we obtain fl3| that typical momentum devia¬ 
tion from a hot spot is ~ 0.06 x 2i r/a. This is a fairly 
small momentum range. As a comparison, for the disper¬ 
sion taken in [l|, the separation between neighboring hot 
spots is much higher: 0.2 x 2ir / a. In this sense already 
for £ = 3a, SC and CDW instabilities come from vicinity 
of hot spots, and this is even more so for £af = 5a, 10a 
also considered in 0 . 

To summarize, we have shown explicitly that anti¬ 
nesting for a half of hot spots does not prevent the in¬ 
stability towards CDW-x order (the one with momentum 
(Q, 0) or (0, Q )) as the corresponding eigenvalue X c d w still 
diverges logarithmically at low T. We believe that such 
an instability was not found in the numerical analysis 
in Ref. 0 because (i) these authors used a much larger 
value of Landau damping than the one obtained in self- 
consistent calculations, and (ii) they neglected fermionic 
self-energy and as a result overestimated the reduction of 
T c dw by FS curvature. 

We thank M. Norman for fruitful discussions. The 
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